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Abstract 

A significant discrepancy exists between experiment and calcula- 
tions of the muon's magnetic moment. We find that standard formu- 
las for the hadronic vacuum polarization term have overlooked pio- 
nic states known to exist. Coulomb binding alone guarantees 7r+ 
states that quantum mechanically mix with the p meson. A simple 
2-state mixing model explains the magnetic moment discrepancy for 
a mixing angle of order a ~ lO^'^. The relevant physical state is pre- 
dicted to give a tiny observable bump in the ratio K( s ) of e^e^ an- 
nihilation at a low energy not previously searched. The burden of 
proof is reversed for claims that conventional physics cannot explain 
the muon's anomalous moment. 

1. Calculations of the muon's magnetic moment do not currently agree 
with experiment. The discrepancy is of order three standard deviations 
and quite important. Among other things, uncertainties of the anomalous 
moment feed directly to precision tests of the Standard Model, including 
the Higgs mass, as well as providing primary constraints on new physics 
such as supersymmetry. In terms of the Lande g factor, the "anomaly" 
a^i = {g — 2)/2 experimentally observed in muons has become the quintessen- 
tial precision test of quantum electrodynamics ( QED ). The current world 
average for fl^ is 11113 : 

^expenmental ^ ( n659208.0 ± 6.3 ) X IQ-'^. 
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The Standard Model theoretical prediction [H for is 

^theory ^ ^ 11659180.4 ±5.1) X IQ-^^. (1) 

The values include contributions from QED, weak interactions including 
W and Z exchanges, plus hadronic ( strongly interacting ) effects. Uncer- 
tainties treat the errors for the hadronic, light-by-light scattering and elec- 
troweak contributions in quadrature. The comparison to 9 digit accuracy 
represents thousands of physicist-years of experimental and theoretical la- 
bor. 

We reviewed contributions considered the most problematic. The hadronic 
vacuum polarization term involves virtual photons emitted by the muon, 
which fluctuate into strongly interacting particles before being re-absorbed. 
The corresponding contribution a^'^' is expressed using the ratio R{s) = 
(r(e+e~ hadrons) / a{e^ e~ ), where a represents the total 

cross section. Because R{s) involves strong interactions it must be taken 
from experiments. The formula for a^'"^' ^'^ is H 



MLO ^ __ I ds "-"-^-^ R( s ), (2) 



^2 p K{x{s)) 

>/s„„-„ S 

where K{x) is a certain QED kernel [4 J given by 

K{x) = + + (^ln(l + x)-x + y 

1 + X 2, 

+ X Inx, 

1 - X 



and x{s) = {1 - / {1 + and = y ^ ~ ^mj/s. We noticed that 

current literature imposes the lower limit of integration Smin = 4m^ in 
Eq. 121. where m„ is the charged pion mass. The threshold comes from as- 
suming the lowest energy hadronic final state is a pion-anti pion pair at 
rest. 

Eq. |2] originates in dispersion relations for complex analytic functions. 
Setting y/Smin = 2m7r ^ 279 MeV is an error which neglects states below 
threshold. The magnitude and sign of the fl^, discrepancy is a shortfall of 
about 28 xlO"^*^ units, relative to some 690 xlO"^*^ units integrating s > 
4m^, suggesting that neglected states might account for the discrepancy. 
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Both continuum and bound states of discrete energy must be taken into 
account. 

Among continuum states there are several contributions, but all ap- 
pear small. Pions are unstable particles for which "anomalous" thresh- 
olds exist, a tiny effect. Next, the production of n'^ plus photon ( 7 ) starts 
at i/s > 135 MeV. The famous "axial anomaly" reasonably predicts the 
7^0 _^ coupling. This channel has been previously noted in the region 
well above threshold, and for 500 MeV < ^/s <1.8 GeV it contributes ||5l 
0.93 10~i° units to aj'''^'^^. Simple estimates find the sub-threshold contin- 
uum "anomalies" too small for our "anomaly. " 

We turn to neglected bound states. The Coulomb interaction will form 
pionic "atoms" bound in the long distance region. The bound states con- 
tributing to vacuum polarization have quantum numbers J^*^ = 
where / is the spin, P the parity and C is charge conjugation. Both isospin 
J = 0, 1 are allowed. A spin-1 "P-wave" pionic state contributing directly 
to vacuum polarization undoubtedly exists. In order to establish a dis- 
crepancy between conventional physics and calculations, one must first 
resolve the contribution of the omitted state( s ). 

It may appear simple to calculate the spin-1 state's contribution to a^i. 
There are several ways to make an estimate. One method approaches the 
problem perturb atively, and extrapolates the Coulomb atomic wave func- 



tion into the core region. Following the treatment outlined in Ref. [|13|, 
the decay width from a P— wave state of pointlike, spinless constituents 
is given by 

\R'( r)]\2 

r( spin - 1 ^ e+e- ) = 8^^ ' \ ^' . (3) 

Here a. is the fine structure constant and R'{0) is the derivative of the 
Schrodinger radial wave function at the origin. There is an important issue 
of what might be meant by the "Schrodinger radial wave function at the 
origin" for pions, which have finite size, but we continue. In the Coulomb 
model R'{0) - ( txifiji ^ ( see below ), and the coupling to 6^6 is far 
too small to matter. For another estimate recall that the p meson mixes 
with the photon at the order of a. The vector, isovector bound state of 
pions must also mix with the p and 7, giving no reason to forbid mixing of 
order 1% or larger, especially given the large pnn coupling. This estimate 
then produces a contribution too large to explain the data for a^i. 
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Why do basic estimates produce such large variation? In Eq. |3] we 
see ( oc'^'^ )^ ^ 10~^^ rules the perturbative Coulomb model, which has 
very little wave function near the origin. Strong interactions and p mix- 
ing, however, produce all their effects in the 1 Fm region of the origin. 
Since the Coulomb model predicts nearly "zero," then neglecting effects 
in the strongly interacting core of the state produces an enormous relative 
error in the contribution to a^i . The only way one might confidently use the 
"Coulomb core" model would be to know in advance that strong interac- 
tion effects and p mixing are even smaller than the electromagnetic binding 
effects, which we find absurd. 

Since strong interactions of hadrons at short distances are a non-perturbative 
problem it would be unwise to proceed on model estimates. The physical, 
strongly interacting mixed state needs a particular name: We call the spin- 
1, isospin-1 state pi-rhonium, symbol n2/p- 

For a model-independent approadi0 we investigated the coupling needed 
for pi-rhonium to be relevant to aj" '^'^ . It is safe to assume the width 
r( 7r2/p e^e~ ) is small and treat its effects as an isolated spike ( delta 
function ) in the integration. A short calculation gives 
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If the state accounts for the full discrepancy indicated by Eq. [T] we then 
have r( 712/ p ) = 28 eV. 

It is interesting that the energy region to detect 7i2/p appears never to 
have been searched. The energies involved seem to have been too low to 
be "interesting." The first e^e~ colliding beams at Frascati [lOJ apparently 
searched well above pion pair threshold to assure data would exist. Data 
available [jllj on e^e~ 7i7i begins at energies well above threshold: 
the high precision CMD-2 experiment lfT2| recently reported on ^/s > 370 
MeV. 

The information can be reversed to predict the effects on a{ e^e~ 
hadrons ). Fig. [T] shows a tiny bump in J^(s) that should be observed at 
energies just below Im-n- The observable effect in i^^'^'^^ is non-negligible 
because the convolution with K{s)/s ( Eq. |2] ) greatly exaggerates the 



^Nothing in this approach hinges on the precise value of binding energy, which needs 
to be found experunentally. 
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Figure 1: Data l|TT| for the ratio R{s) = cr{ e+e^ hadrons) / a{ e+e^ ) 
versus center of mass energy a/s with peaks due to the p, co and (p mesons. The 
contribution of the state 7Z2/p is predicted to be the small bump below thresh- 
old to produce pions. The rcz/p width has been convolved with 10 MeV energy 
resolution. 



small s region. The w^idth of the bump shown includes a 10 MeV figure 
for the energy resolution of the beam and detectors, based on performance 
of modern e^e~ colliders operating in the several GeV region. 

Assuming energy resolution scales in proportion to system energy, the 
calculation indicates that devices operating with current technology in the 
200-300 MeV region can discover n2/p- Signals include the e^e~ peak, 
and a channel we predict with r( n2/p ) ^ 0.6 r( 112/ p 

e^e~ ), based on the phase space. Strong decays will also give signals dis- 
cussed below. Other experiments using rare particle decays, or hadronic 
beams on fixed targets are also attractive: the spin-1 nature of the state 
gives well-established angular distributions to decay products that pro- 
vide distinctive signals. 

This concludes the "IF... THEN" logic of our analysis, whose central 
elements do not depend on any model at all. For example if 35% of the g — 
2 discrepancy comes from hadronic light-by-light scattering, the predicted 
width will be revised. Our analysis also changes the role of "new physics" 
contributions. At most constrains the sum of 7i2/p and new-physics 
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contributions, opening the door for a wider range of new physics than 
constrained to obtain exactly with fine-tuned parameters. ( There is no 
"fine-tuning" of 712/ p properties as the remainder, because the observed 
discrepancy of is defined by what is experimentally measured. ) IF 
the state 7i2/p is observed in e~^e~ channels within error bars of the width 
estimated, THEN the anomalous moment problem is solved. 

2. The model-independent analysis is complemented by exploring rea- 
sonable models in more detail. It is interesting to see what strong interac- 
tion phenomenology may give, while keeping in mind that experimental 
resolution is invariably the only resolution of non-perturbative strong in- 
teraction questions. 

Wave Function Formalism: Suppose one works with Schrodinger wave 
function formalism. Most previous studies of the S-wave pionic atom, 
called "pionium", adopt this approach. 

The literature on pionium is very old, and comes in two varieties. There 
is a long line of research on the bound state physics. We need the effects 
of the short-distance strong interactions on the derivative of wave func- 
tions. Suebka's solution to the S-wave is illustrative[8J. It shows a jump 
of relative order unity in a distance of order l/mn- The average deriva- 
tive Alp/ Ar changes by a factor of about 1/a ^ 10^. This reiterates that 
pi-rhonium's contribution to fl^, is exquisitely sensitive to the model of the 
wave function due to three effects: first, the characteristic size Ax changes 
enormously from Bohr radius-level to Fm-scale level ( a factor of order 
390 ); second, the change of the P — wave is compounded by the kinematic 
dependence to the 5/2 power, and third the wave function is squared to 
get the coupling contributing to a^i . Finally the wave function of a poten- 
tial model cannot be disentangled from how the wave function is used in 
mixing, discussed shortly. 

We could find no discussion of P-wave 1 = 1 pionium in the literature. 
This is perhaps because the mixing problem is hard, and pionium work 
focuses on the scalar state decaying to noTZQ. The P-wave non-relativistic 
model has a centrifugal barrier, but the effective centrifugal barrier of 
strongly interacting systems depends on the potential model. In the Klein- 
Gordon Coulomb problem, which can be solved exactly [T4|, the non-relativistic 
centrifugal barrier l{l + 1) /r^ ^ 1(1 + 1) /r^ + V^{r) + 0{V ) where / is 
the orbital angular momentum and V{r) is the potential of a vector inter- 
action. This is a textbook consequence of gauge invariant vector couplings. 
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Figure 2: Illustration of model sensitivity in a toy model: the ratio of 
space derivatives of Klein-Gordon- Yukawa P-wave functions compared 
to the non-relativistic Coulomb P-wave function. Couplings gpnn equals 
(blue), 3.53 (red), 3.74 (yellow), 4.21 (green), 4.33 (light blue). Fat vertical 
line shows the region r ^ l/niji. 



We found it interesting to explore numerical solutions ( Fig. |2] ) of a toy 
Coulomb plus Yukawa model with coupling ccg = g^rcji/^TC. Range l/niTj 
minimizes parameters and serves as something representing the pion size. 
It is well known| [T4l[T5l[T8| that sufficiently strong potentials overcome the 
non-relativistic centrifugal barrier. Solutions become unstable at coupling 
gpnn > 4.2, oig > 1.4. ( In fact the P-wave ) diverges when literally 
evaluated at zero regardless of coupling strength. ) A large literature on 
relativistic bound states of Klein-Gordon and Dirac type explores many 
variations, including the interesting f act|[18!] that regulating the 1/r singu- 
larities with smeared out densities does not remove strong-coupling insta- 
bility. Softening of the centrifugal barrier is not only a feature of the Klein- 
Gordon system, it is seen in the Bethe-Salpeter formalism[17J and in the 
spinless Salpeter equation[16J that uses Hamiltonian H = ^/p^ + + V. 
One can evidently fit a value of R' ( ) with a defensible range of gpnn- Yet 
the importance of the anomalous moment problem is too high to let nn 
interactions models decide it one way or the other. 

Are there experimental constraints from pionium? Much pionium re- 
search fixes the isosinglet S-wave function to the Coulomb one, and seeks 
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to calculate the pionium lifetime via annihilation tC^tC n^n^. Note 
that mixing for the scalar isoscalar ( S-wave ) state is a much different 
issue than the vector-isovector state we discuss. There is no model inde- 
pendent way to go from the S-wave to the P-wave. Recently the DIRAC 
experiment[6j reported observation of pionium in experiments sending 
particle beams through metal foils. The experiment did not measure the 
spin or isospin of the state. What was actually measured were certain cor- 
relations of longitudinal pion momenta. The correlations were interpreted 
as a lifetime after subtracting a Monte Carlo model of propagation and 
backgrounds. If the P-wave couples substantially, as indicates, the in- 
terpretation of what was observed might be revised. There are too many 
assumptions and too many steps for the limited observations made so far 
to say anything about the 712/ p contribution. 

Mixing Models: Whether or not one might estimate the wave function 
of the spin-1 pionic state on its own, quantum mechanical mixing creates 
new complexity. We explored a simple two state mixing model with off- 
diagonal mixing terms. We mix isovector states and ignore the co meson 
for simplicity. 

The bare basis is denoted \ ft2) , \p), whose linear combinations make 
the physical 1 7r2/p) , \p): 

1^2/ p) \ = ( coS|6 sin^ \ / 1 712) 
\p) J V -s^'"^ cos^ J \ \p) 

The Hamiltonian matrix elements in the bare basis are 

Hint = 



rrift^ mix 
mix mp 



Here mp and mft^ are bare masses, and the mass mixing term is named mi 
Diagonalization gives the mixing angle /3, 



tanjB = mp - mn^ - ^ {mp - ma^ Y + 4m2. ). 

To proceed we must fix three parameters mp, mft^, and mj^. Eigenvalues 
yield the masses of the physical states, namely m^^^^ and mp. Setting the 
physical mass mp = 775.5 MeV fixes one parameter. 

Continuing, we need reference values for the bare states. In quantum 
field theory the calculation will go through a Bethe-Salpeter wave function 
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and irreducible kernels connecting the pionic to photon degrees of free- 
dom. Rather than clutter the calculation with definitions, we express these 
concepts in units of the non-relativistic Schrodinger model with Coulomb 
potential. Notice this is the smallest wave function we might use. No 
particular belief in a Coulomb core is required, because other wave func- 
tions are accomodated by re-scaling. For Hydrogen-like atoms the bind- 
ing energy is given by = —^cp- 1 ( ), where n ^ 2 is the principal 
quantum number in this case. This predicts a reference value for bare mft^^, 
fixing a second parameter. 

The hydrogenic wave function for the 2P state is 

with Bohr radius = 2 / ( (XMji^ and the reduced mass — ntji 
normalization is jQ°°drr-^|R(r )p = 1. The derivative of the bare P- wave 

pionium wave function at the origin is then ^j^^iS^ ) — '^^'^^'^'^■^-'^ • 
consistency of this toy calculation we use the same framework for the 
channels of the p, with width T{ p ^ e+e" ) = 4.70 10"^ Yp = 6.8 keV. 

Finally we use the mixing of wave functions to calculate the width 
r( 7i2/p — > e^e~ ) ^ 28 eV. It fixes the last parameter, so none remain 
free. Since the bare pionium and bare p are orthogonal states, we have 

|J^',^/^(0)|2 = cos2/3|R'^^(0)|2 + sin2/3|4(0)|2; 
|R;(0)|2 = sin2^|R'^^(0)|2 + cos2^|R^(0)|2. 

Solving the equations gives 

sinjS = 0.81%; = 0.46°. 

The wave function of 7t2/p need only be mixed by order the fine structure con- 
stant OL = 1/137, which we find satisfactory. It is just the same order of 
magnitude of mixing one might guess on general grounds. 

With no freedom left the physical mass is mn2/p = 279.11Mey. Several 
digit accuracy is reported for completeness and to exhibit small binding 
effects. The estimated binding energy is about 32 keV, or 65 times the 
"Bohr-atom" estimate. The mixing interaction energy mix — 3-96 MeV. 
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This coupling at order ^ a. x 2m ji seems at the lower limit of what is rea- 
sonable. For one thing, finer discrimination would need control of electro- 
magnetic ( 71+ — ) mass differences. 

These calculations are done in units of the Coulomb core for the bare P- 
wave. Keeping the width of 712/ p fixed by data, larger core wave functions 
will decrease the mixing angle correspondingly. It should be clear the idea 
of wave functions of pions in the p is somewhat naive, and that much 
more information about the short-distance nature of the p is needed to 
believe wave function or mixing calculations. The wave function and the 
mixing ambiguities are complementary "microscopes" for looking into the 
interior of the state. In the core we believe that quarks are relevant, and 
they suggest an interesting model summarizing model requirements. 

In the quark basis production of pi-rhonium goes through a pointlike 
photon vertex where quark chirality is conserved for light u, d quarks with 
masses of a few MeV. The quark-antiquark pair when produced have spins 
aligned in an orbital S-wave. That is, the quark wave function literally at 
the origin is just like that of a ,0. This is a sure fact of short-distance QCD 
perturbation theory. Since the quarks are in an S-wave, no centrifugal 
barrier exists in the interior of the 1 Fm confining distance scale. 

If the mass is as predicted, the 712/ p may be lighter than the S-wave, and 
therefore the second lightest hadron, coming right after the pion triplet. 
This causes a minor puzzle. It is sometimes thought that the lowest state of 
any bound system is an S-wave. The precise statement is that the ground 
state must have no nodes. Feynman[19J gives a simple proof that does 
not depend on the Hamiltonian. The overall quark-model wave function of 
space <S) spin is then a mixture of S and P waves and has no nodes. As the 
quarks emerge into the Bohr radius regime, they are dressed into pions, 
converting quark spin to orbital angular momentum of the pions. When 
pions return to short distance they re-convert orbital angular momentum 
to quark spin. Consider the same physical picture in pion basis. Let v be the 
velocity of the pions, and estimate the radius r^in where pions separate to 
develop one unit of orbital angular momentum: 

rmin X mjrV - h. 

If the pions are non-relativistic v << 1, and r„„„ >> fi/mji. However 
this contradicts confinement, from which we know pions form in the re- 
gion r < h/m-rc- This tells us one cannot extrapolate any non-relativistic 
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P-wave model to Fm distances. The consequences of quark spin-orbital ef- 
fects, in fact, are that neither the tii/ p core nor the p meson can sensibly be 
conceived with pions. Yet pions are perfectly consistent far away from the 
core. We can calculate the velocity Vi^ a from Coulomb physics. Then we 
have Tynin h.lvQmTi '^1 {^ocmji), which is just the Bohr radius r^fn ~ aq- 
If the physical picture seems a bit visionary, the p meson accomplishes an 
identical conversion of quark spin to pion orbital angular momentum in 
every strong decay p nn. This is why we know it applies. The strong 
core of TTiip is a tiny effect. It has an overall probability szn^/S < 10"'*. Yet 
fl^ is a special observable that is highly sensitive to small probability due 
to the close coincidence of the muon and pion masses. 

There are also time scale issues to check. Confinement shows that 
quarks convert to pions in a time of order 1 Fm/c, which is fast com- 
pared to the bound state oscillation time scales. The bound state "or- 
bit time" Torhiti is semi-classically estimated from TorUt ^ a^lv^, where 
fl* is the largest spatial size of the state and a characteristic velocity. 
Semi-classical consistency needs Toy\,n < t, where T is the lifetime of the 
state. We cannot calculate the lifetime, but we observe the strong decay 
712/ p ~^ 71*^7 is allowed in the core. We create a rough estimate by taking 
r( 7i2/p — > e^e~ ) = 28 eV and rescaling it by proportions found in the p 
meson. By this Tj^^^^jot 28eVTp^jjj/T(p — > ) x ps, where ps is a 
calculable phase space ratio. Evaluating the formulas gives Torh/T^ decay 
0{ 10~^ ), which is very acceptable. 

The core amplitude of 712/ p, namely the region contributing to the muon's 
anomalous moment, must be a small strongly interacting object that acts 
as a new hadron. Unless deeply bound the mass of 7i2/p scales with ntjj:, 
and vanishes in the chiral limit. It is very interesting to ask if physics of the 
pion as a Goldstone boson predicts any massless vector particles. The non- 
linear sigma model uses a derivative expansion of the couplings of unitary 
fields U. Within that context, massless vector mesons have long existed in 
gauge-invariant low-energy derivative expansion taking the form 

L = 1/2 fr( a^U - ) ( dj,U -Vj, )')- tr{ F^'Fj.y ) + ... (4) 

Here tr indicates the trace, is the vector combination of UdiiU^ and F^v 
is its particular gauge field strength. Traditionally the gauge-like fields are 
interpreted as the p meson. The need to give the p meson a mass then 
forces ad hoc breaking of the symmetries of the theory. Yet we have 7Z2/p a 
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candidate massless meson in the chiral limit. Was it overlooked twice? It 
is not necessary that Ti2/p remain in the spectrum as electromagnetism is 
turned off. The sigma model is so restrictive that a vector state couples as 
Eq. m predicts as electromagnetism is turned on. Certainly 7i2/p will mix 
with the effective vector fields of the chiral theory. And how much does it 
mix? This is a fascinating avenue, with the implications that the couplings 
of soft pions to 712/ p might be theoretically interesting. 

3. To summarize: A light pionic state with ]^'^ = 1~~ is known to exist 
and was long overlooked. Its contribution to is significantly enhanced 
due to kinematic factors in the basic formula of Eq. |2l The state stands to 
explain the muon magnetic moment anomaly if the p mixing amplitude is 
of order 1%. This is a modest and credible requirement, given the uncer- 
tainties of strong interaction phenomenology. We have made more than a 
good faith effort to check reasonableness of the partial width and its role 
in a^'^' ^'^ . Calculations are non-perturbative and model dependent, and 
far more work would be needed to claim the physical state has the width 
needed. Conversely it would be irresponsible to dismiss the contribution 
of 712/ p ori the basis of wave function or mixing models at their current 
level of development. 

The burden of proof is reversed in our view. If there is truly a discrep- 
ancy between QED and experiment, it seems necessary for proponents 
of the discrepancy to resolve the pi-rhonium contribution. The door is 
opened for new physics to explain part of the discrepancy without obli- 
gation to explain all the discrepancy. The most productive resolution will 
come by comparing experiment to experiment. If pi-rhonium accounts for 
the fl^ discrepancy, there are observable consequences in the experimen- 
tally observable couplings of 7i2/p e^e~ , , 7i^^ 111, 
and similar channels. 
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